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*5 ' Abstract 

We study constraints on the models of bosonic super-weakly interacting particle (super- WIMP) dark 

d : 

. matter (DM) with DM masses mx ~ 0(1 — 100) keV from leptonic decays M — > £u£ + X, where 
,_, ! M = B ± ,D ± ,Df is a heavy meson state. We focus on two cases where X denotes either a light 
' pseudoscalar (axion-like), or a light vector state that couples to the standard model (SM) through kinetic 

o : 

mixing. We note that for a small DM mass these decays are separately sensitive to DM couplings to 



quarks, but not its mass. 
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I. INTRODUCTION 



There is evidence that the amount of dark matter (DM) in the Universe by far dominates that 
of the luminous matter. It comes from a variety of cosmological sources such as the rotation curves 
of galaxies [lj], gravitational lensing Q, features of CMB 3] and large scale structures 4]. While 
the presence of DM is firmly established, its basic properties are still subject of a debate. If dark 
matter is comprised of some fundamental particle, experimentally-measured properties, such as 
its relic abundance or production cross-sections can be predicted. Experimental measurements of 
the abundance VL^m^ 2 ~ 0.12 by WMAP collaboration [5] can be used to place constraints on the 
masses and interaction strengths of those DM particles. Indeed, the relation 



M 2 



^DAih 2 ~ (cr ann v re i) 1 oc — -, 



9' 



with M and g being the mass and the interaction strength associated with DM annihilation, 
implies that, for a weakly-interacting massive particle (WIMP) of DM, the mass scale should be 
set around the electroweak scale. Yet, difficulties in understanding of small-scale gravitational 
clustering in numerical simulations with WIMPs may lead to preference being given to much 



lighter DM particles. Particularly there 
particle with masses of the keV range 



las been interest in studying models of light dark matter 
. According to Eq. (pQ), the light mass of dark matter 
particle then implies a superweak interaction between the dark matter and standard model (SM) 
sector Several models with light O(keV-MeV) DM particles, or super- WIMPs, have been 
proposed jf], Q]. 

One of the main features of the super- WIMP models is that DM particles do not need to be 
stable against decays to even lighter SM particles jg]. This implies that one does not need to 
impose an ad-hoc Z 2 symmetry when constructing an effective Lagrangian for DM interactions 
with the standard model fields, so DM particles can be emitted and absorbed by SM particles. Due 
to their extremely small couplings to the SM particles, experimental searches for super- WIMPs 
must be performed at experiments where large statistics is available. In addition, the experiments 
must be able to resolve signals with missing energy 9[. Super-B factories fit this bill perfectly. 



In this paper we focus on bosonic super- WIMP models 



7| for dark matter candidates and 
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attempt to constrain their couplings with the standard model through examining leptonic meson 
decays. The idea is quite straightforward. In the standard model the leptonic decay width of, say, 
a S-meson, i.e. the process B — > £v, is helicity-suppressed by (m^/m^) 2 due to the left-handed 



nature of weak interactions 



io| 



ns -+ iv) = ^ik,i 2 /k4(i - 4V- (2) 

87r m B \ m B J 

Similar formulas are available for charmed meson D + and D s decays with obvious substitution 
of parameters. The only non-perturbative parameter affecting Eq. ([2]), the heavy meson decay 
constant can be reliably estimated on the lattice so the branching ratio for this process 
can be predicted quite reliably. 

The helicity suppression arises from the necessary helicity flip on the outgoing lepton due to an- 
gular momentum conservation as initial state meson is spinless. The suppression can be overcome 



by introducing a third particle to the final state that contributes to total angular momentum [11 1 
(see Fig. [TJ). If that particle is a light DM candidate, helicity suppression is traded for a small 
coupling strength of DM-SM interaction. In this case, the charged lepton spectrum of the 3-body 
B — > £ug + X (with X being the DM candidate) process will be markedly different from the 
spectrum of two-body B — > Ivi decay. Then, the rate for the process B — > I + p, with $ being 
missing energy, can be used to constrain properties of light DM particles. 

We shall consider two examples of super- WIMPs, the "dark photon" spin-1 particle, and a 
spin-0, axion-like state. The discussion of the vector dark matter effects is similar to a calculation 



of the radiative leptonic decay [llj] , i.e. the spin of the added DM particle brings the required unit 
of angular momentum. In the case of axion-like DM candidate, there is a derivative coupling to 
the SM allowing the pseudoscalar particle to carry orbital angular momentum and hence overcome 
helicity suppression as well. As a side note, we add that the models of new physics considered 
here are very different from the models that are usually constrained in the new physics searches 
with leptonic decays of heavy mesons [l2 1 . 

This paper is organized as follows. In Section [Til we examine the decay width for the process 
M — > Ivt) + X for X = a being a spin-0 particle. We consider a particular two-Higgs doublet 
model, taking into account DM-Higgs mixing in Section [TTTI In Section [IV] we consider constraints 
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FIG. 1: Diagrams for the super-WIMP emission in B — > IvgX. Similar diagrams exist for decays. 
Note that the graph (b) is absent for the vector light dark matter particles discussed in section IIV1 

on a spin-1 super-WIMP candidate. We conclude in Section IVl 



II. SIMPLE AXION-LIKE DARK MATTER 



We consider first an "axion-like" dark matter (ALDM) model, as suggested in [6] and study the 
tree-level interactions with the standard model fermions. The most general Lagrangian consists 
of a combination of dimension-five operators, 

C a = -^7 M 75^ + %F, U F^, (3) 

J a J a 

where X = a is the DM particle and the coupling constant f a has units of mass. Taking into 
account the chiral anomaly we can substitute the second term with a combination of vector and 
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axial-vector fermionic currents 



(\ 4ttC 7 \_ j u , . /8vrC 7 \ - 

£ a = - "T + -* — ^ a VTlhV - %m i> —r — 1 «^75^- (4) 
\fa J a a J V f a a J 

The Feynman diagrams that contribute to the meson decay, for example B — > tvn + a, are shown 
by Fig [TJ The amplitude for the emission of a in the transition M — > tvi + a can be written as 

Au-^iva = Ag + A g , (5) 

where A q , the quark contribution, represents emission of a from the quarks that build up the 
meson and At, the leptonic contribution, describe emission of a from the final state leptons. 

Let's consider the lepton amplitude first. Here we can parameterize the axial matrix elements 
contained in the amplitude in terms of the decay constant fg such as 

(0\u^ l5 b\B) = tf B P£, (6) 

If the mass of the axion-like DM particle is small (m a — > 0), the leptonic contribution simplifies 
to 

At = iV2G F V ub ^-m e ( ™ e [u t Jk{\ - 75)^] - [u e {l - 7 5 )^]V (7) 
f a \2k-pe - J 

Here k is the DM momentum. Clearly, this contribution is proportional to the lepton mass and 
can, in principle, be neglected in what follows. The contribution to the decay amplitude from the 
DM emission from the quark current is 

A q = i (0\uT' i b\B)[utJn(l - 7sK] (8) 
where the current uT^b is obtained from the diagrams in Figured] (a) and (c), 

(9) 



G F 

-Vub 



V2fa 



(^75) (A- fan + m u )Y(l - 7s) 7 M (1 - 7s) (A- h + ™-b){fah) 



m\ — 2p u ■ k ml — 2p b ■ k 



Since the meson is a bound state of quarks we must use a model to describe the effective quark 



antiquark distribution. W'c choose to follow Rets. and 



15], where the wave function for a 
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ground state meson M can be written in the form 



ipM = -|0m(^)75(/Pm + M M g M (x)). 



(10) 



Here I c is the identity in color space and x is the momentum fraction carried by one of the quarks. 
For a heavy meson H it would be convenient to assign i as a momentum fraction carried by the 
heavy quark. Also, for a heavy meson, gjj ~ 1, and in the case of a light meson = 0. For the 
distribution amplitudes of a heavy or light meson we use 



i>£ ~ x(l — x) 

,2 /ji#2 ' 



hi 



(m 2 /M 2 H ) 
1 — x 



I-i 



n -2 



'111 



(12) 



where m is the mass of the light quark and the meson decay constant is related to the normalization 
of the distribution amplitude, 

(13) 



1 A I \A f M 

<p M {x)dx 



2y/E' 



The matrix element can then be calculated by integrating over the momentum fraction 15 ] 



(0\qT^Q\M) = f 1 dx Tr[r^ M ] 
Jo 



(14) 



Neglecting the mass of the axion-like DM particle, the decay amplitude in the B ± case simplifies 
to 



Ja{K ■ r B ) ^ 



(15) 



where rrib is the mass of the 6-quark (or, in general, a down-type quark in the decay), and we 
defined 



M 



1 MX \ x n dx 
x(l — X) 



(16) 
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The total decay width is, then, 



^(2p 2 + 3p 4 + i2p 4 logp-6p 6 + p s 
o 



G 2 F fl\V ub \ 2 M% 
64a* ft 

+ g 2 B $(m 6 ,M B ) 2 (l-6p 2 -i2p 4 logp + 3p 4 + 2p 6 )J , (17) 

where p = mi/ tub- Also, 

$(m 6 ,M B ) = — — . (18) 

jb-IVIb 

Note that <&(mb, Mb) oc 1/m, which is consistent with spin-flipping transition in a quark model, 
which would explain why this part of the decay rate is not proportional to mg. Similar results 
for other heavy mesons, like D + and Df are obtained by the obvious substitution of relevant 
parameters, such as masses, decay constants and CKM matrix elements. 

Experimentally, the leptonic decays of heavy mesons are best studied at the e + e~ flavor factories 
where a pair of M + M~ heavy mesons are created. The study is usually done by fully reconstructing 
one of the heavy mesons and then by finding a candidate lepton track of opposite sign to the tagged 
meson. The kinematical constraints on the lepton are then used to identify the decays with missing 
energy as leptonic decay. 

In the future super-B factories, special studies of the lepton spectrum in M — > £+missing energy 
can be done using this technique to constrain the DM parameters from Eq. (ITT)) . The lepton energy 
distributions, which are expected to quite different for the three-body decays B~ — > ai~i>i are 
shown (normalized) in Fig. [2] for each lepton decay process. However, we can put some constraints 
on the DM coupling parameters using the currently available data on M — > ivg. The experimental 
procedure outlined above implies that what is experimentally detected is the combination, 



T cxp (M iv t ) = r SM (M £p t ) + I dE a 

J E<En 



dT(M -»■ atv t ) 



'e<e dE a 
= T SM (M lv t ) [1 + R a {E )) , (19) 

where Eq is the energy cutoff that is specific for each experiment. Equivalently, cutoff in q 2 can 
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Channel 


Experiment 


Standard 


f%R a (E ) 








(Seen) 


(Maximum) 


Model 


E = 100 MeV 


E' = 50 MeV 


E' = 100 MeV 


E' = 300 MeV 


B(B ± -> T^Vr) 


1.7 x 1CT 4 


7.9 x 10~ 5 


1.6 x 10 2 


4.9 x 10~ 5 


1.9 x 10~ 4 


1.9 x 10~ 3 


B(D ± -> i^vp) 


3.8 x 1CT 4 


3.6 x 10~ 4 


3.1 x 10 3 


4.0 x 10~ 3 


1.8 x 10~ 2 


1.7 x 10~ 2 




5.9 x 10~ 3 


5.3 x 10~ 3 


4.6 x 10 2 


2.0 x 10~ 4 


7.8 x 10~ 4 


6.0 x 10~ 3 




5.4 x 1CT 2 


5.1 x 10 2 


6.5 x 10° 


2.1 x 10~ 5 


8.0 x 10~ 5 


6.2 x 10~ 4 


Channel (Unseen) 














B{B ± -> e ± v e ) 


< 1.9 x 10" 6 


8.3 x 10~ 12 


6.6 x 10 7 


4.6 x 10 2 


1.8 x 10 3 


1.6 x 10 4 


B(B±^/i±P„) 


< 1.0 x 10~ 6 


3.5 x 10~ 7 


1.8 x 10 3 


1.1 x 10~ 2 


4.3 x 10~ 2 


3.6 x 10- 1 


B(D ± -> e^) 


< 8.8 x 10~ 6 


8.5 x 10~ 9 


3.1 x 10 6 


1.9 x 10 2 


7.6 x 10 2 


7.1 x 10 3 


6(L> ± -> r ± S T ) 


< 1.2 x 10~ 3 


9.7 x 10~ 4 


1.0 x 10 1 


1.7 x 10" 3 


7.7 x 10~ 3 


6.2 x 10~ 2 




< 1.2 x 10~ 4 


1.2 x 10~ 7 


9.8 x 10 6 


8.6 x 10° 


3.3 x 10 1 


2.6 x 10 2 



TABLE I: Constraints on f a from various decays. The last three columns represent possible soft photon 
pollution of M — > tvi decays for three different values of photon energy cutoff. 

also be used. In the above formula we defined 

*.(*)= = J tE ^-l^ . (20) 

r SM (M -> tut) Je<e dE a 

Our bounds on the DM couplings from different decay modes are reported in Table H for the cutoff 
values of Eq = 100 MeV. Note that similar expressions for the leptonic decays of the light mesons, 
such as 7r — > atv and K — > aiv come out to be proportional to the mass of the final state lepton. 
This is due to the fact that in the light meson decay the term proportional to g vanishes. Thus, 
those decays do not offer the same relative enhancement of the three-body decays due to removal 
of the helicity suppression in the two-body channel. It is interesting to note that the same is 



also true for the heavy mesons i 



a naive Non-Relativistic Constituent Quark Model (NRCQM), 



• 3,01 i 



similar to the one used in Refs. 16, 1171 1 is employed. We checked that a simple replacement 

Pb = Pb, Pu = Pb (21) 

m B m B 

advocated in [16|, [17J is equivalent to use of a symmetric (with respect to the momentum fraction 
carried by the heavy quark) distribution amplitude, which is not true in general. 

Currently, the SM predictions for the B~ — > t~i>t decay for i = fi, e are significantly smaller 
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Quark 


Constituent Mass 


m u 


ooo.o ivie V 


m d 


339.5 MeV 


m s 


486 MeV 


m c 


1550 MeV 


m b 


4730 MeV 



TABLE II: Constituent quark masses [221 ] used in calculations, 
than the available experimental upper bounds 



18 



20] , which is due to the smallness of V u b and the 



helicity suppression of this process. Thus, even in the standard model, there is a possibility that 
some of the processes B~ — > ■y s £~i / e, with 7 S being the soft photon , are missed by the experimental 
detector. Such photons would affect the bounds on the DM couplings reported in Table [D 

The issue of the soft photon "contamination" of B~ — > £~i>£ is non-trivial if model-independent 
estimates of the contributions are required (for the most recent studies, see 19]). In order to take 
those into account, the formula in Eq. f fT9|) should be modified to 



Y^M^lut) = T SM (M^£u e ) [1 + R a (E ) + R ls (E' )] . 



(22) 



In general, the experimental soft photon cutoff E' Q could be different from the DM emission cutoff 
Eq. Since we are only interested in the upper bounds on the DM couplings, this issue is not 
very relevant here, as the amplitudes with soft photons do not interfere with the amplitudes with 
DM emission. Nevertheless, for the purpose of completeness, we evaluated the possible impact of 



undetected soft photons using NRCQM as seen in [16|, [17] . The results are presented in Tabled] 
for different values of cutoff on the photon's energy. We present the NRCQM mass parameters in 
Table [II] with the decay constants calculated in 21J. The relevant plots for D (D s ) decays can be 
obtained upon substitution Mb —> Md(,d s ), /b fD(D s ), and V u b —> V c d( cs )- Note that there is no 
CKM suppression for D s decays. In order to bound /„ we use the experimentally seen transitions 
B — > t&, Dt s \ — > fiu, and D s — > tv. We note that the soft photon "contamination" can be quite 
large, up to 10% of the standard model prediction for the two body decay. The resulting fits on f a 
can be found in Table ILTT1 As one can see, the best constraint comes from the D ± — > p^v^ decay 
where experimental and theoretical branching ratios are in close agreement. 



9 



0.0 0.1 0.2 0.3 0.4 0.5 0.0 0.1 0.2 0.3 0.4 0.5 




(c) Df lv t a 



FIG. 2: Normalized electron (dashed) and muon (solid) energy distributions for the heavy (B ,D ,Df) 
meson decay channels. Here m a = and x = E^/mg. 



Channel 


f a ,MeV 


B{B ± -> T^Dr) 


12 




236 




62 
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TABLE III: Constraint on f a using the various seen decay channels. 
III. AXION-LIKE DARK MATTER IN A TYPE II TWO HIGGS DOUBLET MODELS 



A generic axion-like DM considered in the previous section was an example of a simple aug- 
mentation of the standard model by an axion-like dark matter particle. A somewhat different 
picture can emerge if those particles are embedded in more elaborate beyond the standard model 
(BSM) scenarios. For example, in models of heavy dark matter of the "axion portal"-type 23], 
spontaneous breaking of the Peccei-Quinn (PQ) symmetry leads to an axion-like particle that can 
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mix with the CP-odd Higgs A of a two Higgs Doublet model (2HDM). For the sufficiently small 
values of its mass this state itself can play a role of the light DM particle. The decays under 
consideration can be derived from the B — > £h>A° amplitude. An interesting feature of this model 
is the dependence of the light DM coupling upon the quark mass. This means that the decay rate 
would be dominated by the contributions enhanced by the heavy quark mass. This would also 
mean that the astrophysical constraints on the axion-like DM parameters might not probe all of 
the parameter space in this model. 

In a concrete model 23j, the PQ symmetry U(1)pq is broken by a large vacuum expectation 



value (S) = f a ^> vew of a complex scalar singlet $. As in [24|, we shall work in an interaction 
basis so that the axion state appears in $ as 



$ = f a exp 



i a 



V2fa 



(23) 



and A appears in the Higgs doublets in the form 



f«exp 



i cot / 





-A c 



v d exp 



i tan /3 
V^vew 



(24) 



where we suppress the charged and CP-even Higgses for simplicity and define tan/3 = v u /vd and 
vew — + v d = We choose the operator that communicates PQ charge to the standard 
model to be of the form 1 



C = A$ 2 $„$ d + h.c. 



(25) 



This term contains the mass terms and, upon diagonalizing, the physical states in this basis are 
given by [24 1 



a p = a cos 9 — A sin 9 



A° p = a sin 9 + A cos 9 



(26) 
(27) 



1 This is the case of the so-called Dine-Fischler-Srednicki-Zhitnitsky (DFSZ) axion, although other forms of the 
interaction term with other powers of the scalar field i> are possible [24[ . 
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f a (MeV) 




f a (MeV) 


Channel 


tan j3 = 1 


tan/3 = 5 


tan/3 = 10 


tan /3 = 20 


B(B ± -> T^Vr) 


70 


340 


357 


361 




416 


2874 


3078 


3131 




532 


1380 


1499 


1529 



TABLE IV: Constraint on f a using the observed decays for various tan /3s. 

where tan# = (vEw/fa) sin 2/3. Here a p denotes the "physical" axion-like state. Thus, the ampli- 
tude for B — > tva v can be derived from 



M(B -»■ £z/a p ) = - sin -> + cos# -»■ &/a) 



(2f 



In a type II 2HDM 24. |25|. the relevant Yukawa interactions of the CP-odd Higgs with fermions 
are given by 



C 



A°ff 



ig tan (3 T , „ ig cot /3 _ . 
rridd^dA H — m u uj5uA 



(29) 

2miy 2m^ 
where d = {d, s, b} refers to the down type quarks and u = {u, c, t} refers to the up type quarks. 
The interaction with leptons are the same as above with d — > £ and u — > v. 

In the axion portal scenario the axion mass is predicted to lie within a specific range of 360 < 
m a < 800 MeV to explain the galactic positron excess 23[. Using the quark model introduced in 
the previous section we obtain the decay width 



T(B -> 



G 2 F \V ub \ 2 m 3 B 



256tt 3 {f^ + v 2 EW sm 2 2^ 
cos 2/3 (m u $f + m 6 « - $f )) + 5 [m 6 ($f - $*) + m u $f 



n2 



(30) 



12^ log(x fl ) - + 3X 4 + (p - l) 4 (4(p - 2)p + 1) - 12(p - l) 4 log(l - p) 



Here we defined x a = m a /mB, and p = me/ms- If we assume / a 3> f^iy sin 2/3 we can then provide 
bounds on f a as seen in Table [IVl Just like in the previous section, the results for other decays, 
such as D( s ) —> £p£, can be obtained by the trivial substitution of masses and decay constants. 
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IV. LIGHT VECTOR DARK MATTER 



Another possibility for a super- WIMP particle is a light (keV-range) vector dark matter boson 
(LVDM) coupled to the SM solely through kinetic mixing with the hypercharge field strength 6]. 
This can be done consistently by postulating an additional U(l)y symmetry. The relevant terms 
in the Lagrangian are 

C = -F/"" - \v^ v - + ^V^W + C w , (31) 

where Cy contains terms with, say, the Higgs field which breaks the U{l)y symmetry, k parame- 
terizes the strength of kinetic mixing, and, for simplicity, we directly work with the photon field 
Ap. In this Lagrangian only the photon fields (conventionally) couple to the SM fermion 
currents. 

It is convenient to rotate out the kinetic mixing term in Eq. (I3T]) with field redefinitions 

A->A'- = V, V -> H 1 V 1 . (32) 



K 



The mass m v will now be redefined as m v — > J^^ k2 ■ Also, both A'^ and now couple to the SM 
fermion currents via 

C f = -eQfAtfffjf - J^^V$ffi/> f , (33) 

where Qf is the charge of the interacting fermion thus introducing our new vector boson's coupling 
to the SM fermions. Calculations can be now carried out with the approximate modified charge 
coupling for k 1 , 

; « KC 34 

As we can see, in this case the coupling of the physical photon did not change much compared to 
the original field A^, while the DM field acquired small gauge coupling ne. It is now trivial 
to calculate the process B — > EvVdm, as it can be done similarly to the case of the soft photon 
emission in Sect. El Employing the gauge condition e • k = for the DM fields, the amplitudes 
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Channel 


«- 2 i? y (£;o) 

£0 = 100 MeV 


K 


B{B ± 




8.8 x 10~ 3 


< 11.6 


B(D ± 




5.7 x 10- 1 


< 0.31 


B(Df 




5.4 x 10~ 2 


< 1.49 


B(Df 




1.3 x 10~ 4 


< 20.8 


B{B ± 




1.8 x 10 3 


< 11.2 


B(B ± 


-> M ± ^) 


1.0 x lO" 1 


< 4.17 


B(D ± 


-> e ± i/ e ) 


1.5 x 10 3 


< 0.83 


B(D ± 


-> r±P r ) 


1.8 x 10~ 4 


< 36.4 


B{Df 




5.2 x 10 2 


< 1.37 



TABLE V: Constraints on k using various decay channels. All other values are the same as in Table HI 
become in the limit my — > 

'A^{1 - 75 )z/ £ + Bi la {l - l5 )v t + Cj(l - 75 K +D»£<r IM (l + 7s)^l(35) 



_. G F V ub Kee* a 
6k -p B 



with the coefficients 



AH 



B 



3V2f B - 2V / 3(<^ 3 + $f )1 k«q Q - 2V3($£ - 3^f)ie^ ap k a q p , 



3V2f B - 2v / 3($ B + $f )] (A; - g) - -^/BTn| 
2\Z3gm B [m 2 (0 o - 3</>i) + 2771501] , 



3\/2f B m e 
-3V2if B m £ 



k ■ 
k ■ q 



k-p e 



k\ 



(36) 

(37) 
(38) 

(39) 



and q = pt + p u . Again, we fit the parameter k using the same data as in the axion-like DM 
case. The results are shown in Figure fVl where the D — > [l P^V decay can yield the best bound. 
Using the best constraint on k from the D ± — > p^v^V decay we can limit the contribution to 
yet-to-be-seen decays in Table fVTl 

As we can see, the constraints on the kinetic mixing parameter k are not very strong, but could 
be improved in the next round of experiments at super-flavor factories. 
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X 

(c) Df ->■ IviVdm 

FIG. 3: Normalized electron (dashed) and muon (solid) energy distributions for the heavy { (B^ , , Djr } 
meson decay channels. Here m a = and x = Ei/tub- 



Channel 


B(« = 0.31) 


B{B ± -> e*^) 
B(D ± -> e ± P e ) 
B(£>± ^e±P e ) 


1.4 x 1(T 9 

3.6 x 1(T 9 
1.2 x 10" 6 

1.7 x 1(T 8 
6.2 x 1(T 6 



TABLE VI: Contributions to various yet-to-be-seen channels using the the fit on k in Table [V] 
V. CONCLUSIONS 

We considered constraints on the parameters of different types of bosonic super- WIMP dark 
matter from leptonic decays of heavy mesons. The main idea rests with the fact that in the 
standard model the two-body leptonic decay width of a heavy meson M = {B ± , D ± , Df}, or 



15 



T(M — > iu), is helicity-suppressed by (m^/m^) 2 due to the left-handed nature of weak interac- 
tions [10]. A similar three-body decay M — > ivgX decay, which has similar experimental signature, 
is not helicity suppressed. We put constraints on the couplings of such DM particles to quarks. 
We note that the models of new physics considered here are very different from the models that 
are usually constrained in the new physics searches with leptonic decays of heavy mesons [12[| . 

We would like to thank Andrew Blechman, Gil Paz, and Rob Harr for useful discussions. This 
work was supported in part by the U.S. National Science Foundation CAREER Award PHY- 
0547794, and by the U.S. Department of Energy under Contract DE-FG02-96ER41005. 
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